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ABSTRACT 



The Wisconsin H-Alpha Mapper (WHAM) has surveyed the distribution and 
kinematics of ionized gas in the Galaxy above declination —30°. The WHAM 
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Northern Sky Survey (WHAM-NSS) has an angular resolution of one degree and 
provides the first absolutely-calibrated, kincmatically-resolved map of the Ha 
emission from the Warm Ionized Medium (WIM) within ~ ±100 km s~^ of the 
Local Standard of Rest. Leveraging WHAM's 12 km s~^ spectral resolution, we 
have modeled and removed atmospheric emission and zodiacal absorption fea- 
tures from each of the 37,565 spectra. The resulting Ha profiles reveal ionized 
gas detected in nearly every direction on the sky with a sensitivity of 0.15 R (3cr). 
Complex distributions of ionized gas are revealed in the nearby spiral arms up to 
1-2 kpc away from the Galactic plane. Toward the inner Galaxy, the WHAM- 
NSS provides information about the WIM out to the tangent point down to a 
few degrees from the plane. Ionized gas is also detected toward many interme- 
diate velocity clouds at high latitudes. Several new H II regions are revealed 
around early B-stars and evolved stellar cores (sdB/0). This work presents the 
details of the instrument, the survey, and the data reduction techniques. The 
WHAM-NSS is also presented and analyzed for its gross properties. Finally, some 
general conclusions are presented about the nature of the WIM as revealed by 
the WHAM-NSS. 

Subject headings: ISM: structure — Galaxy: halo — HII regions — ISM: atoms 

1. INTRODUCTION 

Studies of faint emission hues over the past two decades have firmly estabhshed the 
Warm Ionized Medium (WIM) as a significant component of the ISM, especially in the 
halos of disk galaxies (Reynolds 1993; Hoopes et al. 1999; Rand 1996; Rossa & Dettmar 
2000; Dettmar 1998). In our own Galaxy, measurements show the WIM to have a mass 
surface density about one-third that of H I, an ionization power requirement equal to the 
kinetic energy from Galactic supernovac, and a characteristic scale height of about 1 kpc 
in the solar neighborhood (Reynolds 1993; Ferriere 2001). Although originally discovered in 
the 1960s through radio observations (Hoyle & EUis 1963; Guelin 1974), information about 
the distribution, kinematics, and other properties of the WIM has been obtained primarily 
through the study of optical emission lines, particularly Ha, [S II] A6716, and [N II] A6583. 
Gharacterizing the details of the WIM is important not only for improving the picture of 
the ISM and the interaction among its components, but also for understanding the nature of 
the foreground gas through which all extragalactic and cosmological observations are made 
{e.g. Lagache et al. 2000; Bennett et al. 2003). 

Although the details of how the WIM is ionized and heated are not well understood 
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yet, many authors have presented ideas for the source of the ionization and heating of the 
WIM and its relationship to other phases of the ISM. Miller & Cox (1993), Dove & ShuU 
(1994), and Dove et al. (2000) present a picture of a pervasive, warm, (originally) neutral 
medium being ionized by radiation from early-type stars leaking between cold H I clouds 
or out of density bounded H II regions. Potentially, such a picture results in regions of 
fully ionized gas and "shadowed" neutral gas. Others suggest that the Ha emission should 
be positively correlated with H I features. Norman (1991) and Koo et al. (1992) associate 
the pervasive Ha emission with the walls of superbubbles, "chimneys," and "worms," while 
McKee & Ostriker (1977) suggest that prominent Ha emission is confined to the surface of 
neutral clouds bathed in a hot medium. Spitzer & Fitzpatrick (1993) and Sciama (1990) 
have suggested that the H"*" should be well-mixed with the neutral material, although these 
scenarios are inconsistent with recent measurements of the ionization fraction within the 
Ha emitting gas (Reynolds et al. 1998; Hansen et al. 2002a). Slavin et al. (2000) have 
predicted a minimum ionizing flux present in the halo of the Milky Way contributed by 
old, cooling supernova remnants. Although this flux is not able by itself to account for the 
ionization of the WIM, they suggest that it may be comparable or even dominant to the 
Lyman continuum flux escaping from the disk at large distances from the plane. In addition 
to contributions from multiple ionization sources, substantial evidence has been provided 
recently that the heating of the WIM is likely dominated at low densities by a source other 
than photoionization (Haffner et al. 1999; Reynolds et al. 1999). 

The development of high quantum efficiency, low read noise CCDs has recently allowed 
breakthrough Ha studies of the WIM. In addition to the velocity-resolved, all-sky survey 
described here, three large-area imaging surveys are also available or in progress. A southern 
survey with arc-minute resolution has been completed by Gaustad et al. (2001). Dennison 
et al. (1998) have been imaging the northern sky with similar resolution. A high- resolution 
(~ 1") survey of the Galactic plane and Magellanic clouds is also underway by Parker & 
Phillipps (1998). Furthermore, direct measurement of extremely faint, optical diagnostic 
hues in the WIM, including He I A5876 (Tufte 1997), [O I] A6300 (Reynolds et al. 1998; 
Hausen et al. 2002a), and [N II] A5755 (Reynolds et al. 2001b), have placed important 
constraints on the nature and source of the ionization and heating of this gas in our Galaxy. 

Some of the first, most comprehensive work on the WIM (also referred to as Diffuse 
Ionized Gas — DIG — in this context) of an edge-on spiral galaxy similar to the Milky Way 
was undertaken by Rand (1997, 1998) who measured numerous faint ionized emission lines far 
into the halo of NGC 891. Domgorgen & Dettmar (1997) and GoUa et al. (1996) also probed 
similar lines in the edge-on galaxies NGC 2188 and NGC 4631. All found a pronounced 
increase in the intensity ratios [S II] /Ha and [N II] /Ha in regions with fainter Ha emission. 
For spectra obtained perpendicular to the galaxies, these ratios rise systematically with 
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distance from the galactic plane. More recent spectroscopic work has extended the detection 

of the WIM and these line ratios to impressive distances above galactic disks (Tiillmann & 
Dettmar 2000; Rand 2000; Colhns & Rand 2001; Otte et al. 2001, 2002; Hoopes & Walterbos 
2003). Meanwhile, deep imaging of the extent and morphology of the WIM in spiral galaxies 
is providing valuable information about the link between the ionized layer and the star 
formation within galaxies (Hoopes et al. 1999; Rossa & Dettmar 2000; Zurita et al. 2000; 
Rossa 2001; Thilker et al. 2002; Zurita et al. 2002). 

Utilizing this new CCD technology, we designed and built a dedicated survey instrument, 
the Wisconsin H- Alpha Mapper (WHAM), which has produced the first deep, velocity- 
resolved maps of the WIM in our Galaxy as well as the first clear detections of Ha emission 
from high- and intermediate-velocity clouds (Tufte et al. 1998; Haffner et al. 1999) and 
detections of extremely weak diagnostic lines from the disk and halo gas of the Milky Way 
(Tufte 1997; Reynolds et al. 1998, 2001b). The WHAM Northern Sky Survey (WHAM-NSS), 
the subject of this work, provides the first Ha survey of the distribution and kinematics of 
the diffuse H II comparable to earher 21 cm surveys of the H I. 

2. OBSERVATIONS 

The WHAM instrument was installed at Kitt Peak in late 1996 and began obtaining 
survey data in early 1997. Much of the survey was complete after this first year. FiUing the 
remaining holes, obtaining backup data for regions that were taken on nights of questionable 
quality, and extending the southern declination limit from —20° to —30° continued through- 
out 1998. Below, we describe the unique features of our instrument (§2.1), the characteristics 
and strategies of the survey itself (§2.2), and the methods we use to transform the raw data 
into clean Galactic profiles (§3). 

2.1. The WHAM Instrument 

Due to the extended nature of diffuse ionized gas, large-aperture Fabry-Perot detection 
techniques have proven to be quite successful in studying the WIM (Reynolds et al. 1990; 
Reynolds 1990). WHAM continues to exploit this feature with its 15-cm aperture, dual- 
etalon Fabry-Perot spectrometer that delivers high-efficiency observations at high spectral 
resolution. WHAM is simply designed and consists of the following primary components, 
which are highhghted in Figure 1: 

• Siderostat. A custom-built all-sky siderostat houses two flat folding mirrors in a 
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tipped "alt-az" mount. The system was designed to point and track continuously over 
the entire sky. The siderostat was designed and built by the University of Wisconsin- 
Madison Physical Sciences Laboratory. 

• Pre-etalon Optics. A pair of 860-cm focal length lenses provides telecentric illumi- 
nation of the etalons. The stationary housing of the siderostat hosts a 63-cm diameter 
lens while a companion 18-cm lens sits in front of the first etalon surface. With this 

configuration, a one-degree beam on the sky is imaged into a 15-cm diameter region 
at the focal plane between the two etalons, and the 0?50 half-angle cone from the sky 
is transformed into a 2?09 half angle cone as it passes through the etalons. 

• Etalons. Two 15-cm etalons form the spectroscopic element for WHAM. The dual- 
etalon design combines a high-resolution etalon (with a gap spacing of 0.0471 cm) 
in series with a low-resolution etalon (with a gap spacing of 0.0201 cm). The low- 
resolution etalon suppresses higher orders from the high-resolution etalon, opening up 
a much larger useable spectral range in each observation than a single high-resolution 
etalon would provide, while suppressing the broad Lorentzian-likc wings of the Fabry- 
Perot response profile. This feature is illustrated in Figure 2. With our chosen gap 
spacings, a reflective finesse of 30, and a maximum half angle of 2?09, WHAM spectral 
observations deliver 12 km resolution (R 25, 000) over a 200 km spectral 
window near Ha. The location of the spectral window is selected by changing the 
pressure of a high index of refraction gas (SFg) in sealed etalon chambers. This change 
in the index of refraction translates the transmission functions shown in Figure 2 
horizontally (i.e. in wavelength). One transmission peak from each etalon must also be 
coincident for maximum efficiency. We perform this alignment nightly as a part of our 
set-up procedure, and the pair of pressures that maximize the composite transmission 
function is called the "tune" for that spectral window. Since environmental effects can 
change the index of refraction as well, the tune can change slowly over day-long spans 
and must be confirmed every night. 

• Post-etalon Optics. To achieve even higher spectral isolation and to reject light 
outside the 4800 A to 8400 A band where the etalons are no longer reflective, the 
final rejection element is a standard 3.1" (7.87 cm) narrow- band (~ 20 A FWHM) 
interference filter with a transmission waveband centered near the emission line of 
interest (see Figure 2). WHAM is fitted with an eight-position filter wheel to allow 
observations of many emission lines during a night. A pair of 104-cm focal length 
lenses provides telecentric imaging of the Fabry-Perot ring pattern that emerges from 
the last etalon surface onto the interference filter. Finally, a four-lens stack (effectively 
f/0.9) images this ring spectrum onto a ~ 1 cm diameter area of the CCD chip, while 
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minimizing aberrations to less than 100//m (one 4x4 binned pixel). 

• CCD. WHAM's imaging device is a cryogenically cooled (~ —101° C) Tektronix 
TK1024 with 24/xm pixels and a quantum efheiency of 88% at Ha. Our detector setup 
results in a gain of 0.57e~/ADU and an effective read noise of 3.87e~ rms after 4x4 
on-chip binning. 

Coatings on the etalons, optics, and CCD allow observations anywhere between 4800 A 
and 7320 A. All observation taken for the survey are centered around Ha and were obtained 
through a 20 A FWHM interference filter manufactured by Barr and Associates whose center 
wavelength is centered near 6565 A on axis, which translates to 6563 A (Ha) at the average 
angle of incidence onto the filter. 

Note that in the configuration described above, WHAM images only the Fabry-Perot 
ring pattern emerging from the etalons. All spatial information is lost, with each point 
within the l°-diametcr beam sampled equally by every spectral element within the 4.4 A 
(200 km s^^) bandwidth of the spectrum. In this "spectral mode," a CCD image contains 
the average spectrum within the one-degree diameter circular beam, avoiding any confusion 
between spectral features and angular structures of the source (or stars) within the beam. 
WHAM also has a set of optics that can be placed in the post-etalon beam to provide 
an extremely narrow-band (20-200 km s"^), 3' angular resolution image of the sky within 
the one-degree field of view. However, only the spectral mode (with one-degree angular 
resolution) is used for the Ha Sky Survey observations presented here. Further details on 
the WHAM spectrometer and its calibration can be found in Tufte (1997). 

WHAM was designed from its inception as a remotely operable facihty. At the core of 
the operational software, WHAM implements the WIYN Telescope Control System (Percival 
1995). All essential aspects of the instrument arc able to be efficiently controlled by a com- 
puter over a low-bandwidth connection to the Internet. With this setup, WHAM observers 
have been able to maintain a demanding observing schedule without the associated costs in 
time and funding of long-distance travel. 

2.2. Survey Strategies 

The northern sky survey consists of more than 37,000 pointings of the one-degree beam 
on a grid with Ab = 0?85, A/ = 0?98/ cos 6, and 6 > —30°. The origin of the longitude grid- 
ding is offset by (0?98/2)/ cos 6 between adjacent latitude rows. This odd spacing achieves 
nearly full sky coverage near the Galactic plane by nesting the beams. At higher latitudes. 
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some "slipping" of the rows begins to occur and at some latitudes slightly larger gaps of 
unsamplcd sky are left between adjacent latitude rows. Since the goal was to probe the 
extended emission of the WIM, these small gaps have little effect on the completeness of the 
survey. To organize these observations, the individual pointings are grouped into observa- 
tional "blocks." Each block consists of up to 49 one-degree pointings, typically in a 7 x 7 
nested grid. An example of our beam spacing and block division is presented in Figure 3. 
Observations begin in one corner of a block and proceed in a boustrophedonic raster toward 
the opposite side, assuring that the sequence of observations within a block takes small steps 
in both time and on the sky. The exposure time for each one-degree pointing is 30 seconds, 
providing a sensitivity level that is better than 0.15 R (3(7). These approximately 30-minute 
blocks then form a tractable observational unit. 

To minimize solar features in the background of our spectra, particularly the Praunhofer 
Ha absorption line, survey observations were obtained only after astronomical twilight had 
ended and when the moon was below the horizon. In addition, all observations are made 
below a zenith distance, ^, less than 60° (i.e. at an airmass, Z, less than 2) to avoid 
substantial atmospheric extinction corrections (see §3.5) and reduce the contamination of 
atmospheric lines (see §3.3.3). 

One final constraint attempts to maximize the spectral velocity separation between the 
geocoronal Ha emission line (the primary atmospheric contaminant in these spectra, see 
§3.3.2) and the Local Standard of Rest (LSR). Due primarily to the orbital velocity of the 
earth and the peculiar velocity of the sun with respect to the LSR, this separation depends 
on both the direction in the sky and the time of the year. Since the bulk of the Galactic Ha 
emission occurs near km s~^ LSR, we chose to observe blocks when the velocity separation 
between the LSR and geocentric frames was > 20 km on a given night. With 12 km s~^ 
velocity resolution, these separations are sufficient to identify reliably the relatively bright 
but narrow geocoronal Ha emission line and subtract it from the spectrum. This criterion 
is satisfied for every direction on some night during the year except near the ecliptic poles, 
where the separation is limited to a maximum of 15 km s~^. 

For most directions in the sky, there are two periods during the year that satisfy this 
constraint, one with the LSR blueward of geocentric zero as the earth's orbital vector points 
toward the general region of the look direction, and one a few months later with the LSR 
shifted to the red as the orbital vector points away from the look direction. Blocks observed 
with the former configuration (LSR at negative geocentric velocity) rise above the horizon 
later in the night, and we label these "morning" blocks. Those observed with the LSR at 
positive geocentric velocities are visible earlier and are labeled "evening" blocks. Although 
WHAM is continuously tunable and can locate the 200 km s~^ window anywhere within the 



-8- 



filter bandpass, we chose to adopt two fixed tunes for these configurations, called the morning 
and evening tunes. Since most of our spectral contamination is atmospheric (and thus fixed 
with respect to the geocentric scale), this decision minimizes variations in the calibration 
phases of data reduction. One negative side effect is that the actual range of LSR velocities 
observed is not fixed from pointing to pointing. As a result, the spectral coverage of an 
individual observation ranges between LSR velocities of [—138, +64] km s^^ and [—71, +131] 
km s""^, with most observations taken within ±10 km s"-*^ of [—116, +87] km s"-*^. The entire 
WHAM-NSS does have complete spectral coverage between —70 and +65 km s~^, however. 
90% of the survey pointings provide complete coverage from —90 to +75 km s~^. 



3. DATA REDUCTION 

Due to the simplicity of the instrument and the relatively stable skies at Kitt Peak 
near Ha, reduction of our spectral data was straightforward. The creation of calibrated 
interstellar Ha spectra from raw WHAM images consisted of the following steps: 



3.1. Ring- Summing 

Figure 4a shows the raw image output from the WHAM CCD camera. This "ring 
image" represents the spectrum within the one-degree beam. To reiterate, unlike many other 
astronomical Fabry-Perot instruments, the spectral mode of WHAM records only spectral 
information — each part of the sky within the beam is sampled equally at all wavelengths 
and therefore no spatial information is obtained. Rings in these images are emission lines 
and the transition to zero signal at the largest radius represents the edge of our circular 
aperture defined by the 0.6-m diameter siderostat lens. Wavelength decreases radially from 
the center to the edge of the ring pattern, which is azimuthally symmetric. Since we are 
imaging only small angles {6 < 2?09) emerging from the etalons, AA ~ f ~ ~ r^, where 
r is the distance of a pixel from the center of the ring image (Tufte 1997). Thus, equal-area 
annuli correspond to equal spectral intervals. 

The image is converted into a one-dimensional spectrum by averaging the signal recorded 
in successive annuli ("annular-summing"; Coaklcy et al. 1996; Nossal et al. 1997). For these 
data, we define an 85-pixel (area) annular window for each data point and average the signal 
of all pixels whose centers fall within the radial hmits for the annulus. This annular window 
size was chosen because it corresponds to a ~ 2 km s~^ spectral window, oversampling the 
width of our instrumental resolution by a factor of about six. The actual velocity width 



-9- 



assigned to each window (equivalently an output spectral data point) has been determined 
empirically through observations of Hydrogen- Deuterium and Thorium-Argon calibration 
lamps. We measured the separation between several pairs of emission lines at different 
positions in the spectral field, allowing us to determine the first and second order coefficients 
of the velocity dispersion. The functional solution for our velocity calibration is then simply 
a second order polynomial: 

Vi [km s"^] =vo + 2.083 v - 0.00036 i^, (1) 

where Vq is an arbitrary zero point (at this stage, see §3.4 below for conversion to the 
LSR frame), v is the integer index of a annular binning window, and Vi is the velocity 
corresponding to that bin. v = corresponds to the outer-most (bluest) annulus. 

Due to the discrete nature of the CCD pixels, these windows will not all have exactly 85 
pixels, but the deviation in bin size is small so that the noise characteristics of the resulting 
data points are similar. We employ 4x4 on-chip binning to reduce read noise, but do 
not compromise spectral resolution at the blue end of the spectrum. The optical setup is 
designed so that the last 8 km s~^ spectral interval is one 4x4 "super-pixel" wide (96 //m 
on the chip) and fully sampled azimuthally. 

Three other minor corrections also occur at this stage: 

• Bias subtraction. Bias frames are taken each night during observations. However, 
since we saw no significant structure in the bias images or drift in the average bias 
value during the first three years of operation, we reduced this step to subtracting a 
bias constant from each image. 

• Cosmic ray removal. Cosmic rays are removed in two steps. The first is a simple, > 
3(7, neighbor-rejection technique that is used in many standard CCD image processing 
packages. We then take advantage of our azimuthal symmetry in the ring-summing to 
do a finer rejection. At a given radius, pixels arc exposed to the same source intensity 
since they lie within the same spectral element. Wc fit a first-order polynomial to the 
/ vs. distribution within this annulus and apply another threc-sigma rejection to 
pixels about this line. In a typical 30 s spectrum, less than 1% of the image pixels are 
rejected with this multi-pass method. 

• Reflection subtraction. A reflection with an intensity of 4% that of the primary ring 
pattern is offset by a; = +37.50, y = —21.25 pixels from the center of the image (very 
faintly visible in the lower right quadrant of Figure 4a). This feature arises from the 
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reflection of the image off a wedged (0.01 rad), uncoated surface in the second etalon. 
To correct for this effect, the raw ring pattern is scaled, translated, and subtracted 
from itself before the ring-summing takes place. Other reflections in the system are at 
a significantly lower level and have no detectable affect on the resulting spectra. 

Figure 4b shows the ring-summed spectrum from a typical WHAM observation. 

3.2. Flat Fielding 

The majority of transmission inhomogcncitics in the WHAM system can be removed by 
taking the spectral ring image produced by shining a diffuse, white light soTircc through the 
entire optical path, processing it into a traditional spectrum as describe above, and dividing 
every spectrum by this resulting "ffat-field spectrum" . Although functionally quite different 
from dividing a spectral image by a white-light image and then ring-summing, we have found 
the former procedure to be much more robust for our application. Since roughly 85 pixels 
contribute to a single spectral data point, we find that the effect of intrinsic pixel to pixel 
gain variations is not significant for our data, especially with our second cosmic ray step. 
Pixels with significantly abnormal sensitivities are automatically rejected with this method, 
improving the estimate of the original source intensity. 

This process effectively removes the primary vignetting shape produced by the geometry 
of the tandem Fabry-Perot etalons and the shift of the transmission function of the narrow- 
band filter with incident angle. Together, these factors result in a ~ 20% reduction in 
transmission from the center to the edge of the ring image. 

3.3. Extracting Galactic Emission 

Three components of our spectra, the background continuum, the geocoronal line, and 
fainter atmospheric lines, are removed in a single fitting and subtraction procedure. They 
are each described in more detail below. Their characterization and subtractions represents 
the bulk of the work to isolate pure Galactic emission in WHAM observations. Custom 
fitting software was designed to be semi-automated during the reduction of an observational 
block so that these corrections would be tractable with our large number of spectra. As 
noted in §2.2, pointings within a block are taken in small steps (~ 1°) across the sky and in 
time (~ 30 seconds) so that reducing many of the blocks usually only requires hand-fitting 
the first pointing of a block. This solution is used as the initial input to the fit for the 
next pointing and the fitting steps along through the block. Galactic emission is modeled 



-li- 



as a series of Gaussian components for this procedure, but we did not focus on the exact 
component structure of the emission in reducing the data. Instead we simply required a good 
fit to the underlying Galactic emission, especially in regions relatively free from atmospheric 
contamination. 



3.3.1. Background Subtraction 

The brightness of the sky continuum emission at Kitt Peak within our one-degree beam 
at Ha is typically 15-30 ADUs, corresponding to approximately 0.022-0.044 R (km s~^)~^ 
(or 1.0-2.0 R A~^), and is extremely flat within our 200 km s~^ bandpass (1 R = 10^/47r 
photons cm" -2 s"^ sterol = 2.4 x 10~^ ergs cm ^ s ^ stcr ^ at Ha). Stellar sources in the 
onc-dcgrcc beam with my < 5 also noticeably elevate the background level in the spectrum. 
A zero- or first-order polynomial is used in the fitting procedure to remove the background. 

If a bright star also has prominent absorption features in the spectral window (typically 
Ha), it can severely contaminate the spectrum and the blending of these stellar features 
with the emission makes recovery of the Galactic emission nearly impossible. In the publicly 
released data set, pointings that contain stars with my < 5.5 (~ 10%) are flagged to denote 
possible contamination. In the maps presented here, the relevant intensity integration of the 
pointing has been replaced by an average from neighboring spectra. WHAM may re-observe 
these directions with a small occulting disk near the imaging focal plane to eliminate the 
stellar contamination, but such observations are not included here. 

3.3.2. Geocoronal Subtraction 

The primary contamination in the survey data arises from the geocoronal Ha emission 
line (see Figure 4b) . This emission is produced primarily from excitation of neutral hydrogen 
in the outer layer of the earth's atmosphere through the scattering of solar Ly/3 radiation 
(n = 1 — > 3). The subsequent decay cascade produces a prominent Ha emission feature 
in all spectra with an intensity of 2-13 R (Nossal et al. 2001). Since the intensity of the 
exciting Ljf3 radiation depends most sensitively on the height of the earth's shadow along 
the look direction, the resulting strength of the contaminating Ha line depends on both the 
direction and time of the observation. For example, at a constant hour angle and declination, 
the intensity of this line can change by a factor of three from a minimum at midnight to 
a maximum at twilight. At our spectral resolution, we do not resolve the fine structure of 
this hue (Nossal et al. 1998) and can fit it well with a single Gaussian component with an 
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intrinsic width of 5-8 km s 



3.3.3. Fainter Atmospheric Line Subtraction 

In addition to the geocorona, several much fainter (< 0.2 R), atmospheric hnes of 
unknown origin are present near Ha in our spectra. As expected, the intensity of such hnes 
rises with increasing zenith distance (i.e. airmass). The hnes appear to be fairly stable over 
at least half-hour periods, but variations from night to night can be as large as a factor of 
two. These lines have been fully characterized by Hansen et al. (2002b), and their measured 
parameters are reproduced here in Table 1. 

To remove the atmospheric emission from the survey spectrum, we created a synthetic 
template of the lines using the positions, widths, and relative intensities listed in Table 1. 

Wc hoped that this template could then be scaled by one parameter to fix the absolute 
intensity of all these lines together. However, an additional parameter was needed to adjust 
the intensity of line 10 at Vgco = +73 km s^^ independently from the rest of the template. 
This "line" is affected by at least two known complications, which likely dictate the need for 
a second parameter in our template. First, the bright OH emission line located at 6553. 6A 
is not fully rejected by our etalon configuration and a faint "ghost" line appears at this 
position near Ha. There is no a priori reason to suspect that the intensity of this ghost 
line is linked to the fainter lines seen near Ha. Second, there is a known absorption feature 
(probably H2O) present at Vgeo = +65 km s~^. The strength of this line is dependent on the 
brightness of the sky background on a given night toward a given direction, and will also 
depend on the column of H2O along the line of sight. Its presence reduces the contribution 
of the other lines near these wavelengths, most substantially, the line at fgco = +75 km s~^. 
In fact, wc find a good correspondence between low values of the atmospheric parameter 
fit to the Vgeo = +75 km s~^ component and higher temperatures outside (when the water 
vapor column in the air can be higher). 

With this template in hand, we created very high signal-to-noisc "block-averaged" spec- 
tra for each survey block (see §2.2) and fitted all the features (Galactic, geocorona, and 
atmospheric) in these spectra to empirically derive the two scaling parameters for the atmo- 
spheric template. Each of these block-averaged spectra represents a spatial average of about 
30-50 square degrees on the sky and a temporal average over about 30 minutes. As noted 
above, the intensity of these faint atmospheric lines changes slowly over hour-long periods 
and in practice, we find that this method removes the contamination from each individual 
survey spectrum quite well. The template of weak atmospheric lines, scaled by these deter- 
mined parameters, is used when determining the geocoronal fit (§3.3.2) and is removed from 
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the Galactic emission profile along with the geocoronal Gaussian. 



3.4. Velocity Registration 



As compensation for the contamination by the geocoronal line, we gain an accurate zero- 
point wavelength calibration in each Ha spectrum. Since the geocoronal line arises primarily 
from ground state excitation to n = 3 and not from recombination, the full assembly of fine 
structure Ha lines are not present and the line center of the geocoronal Ha multiplet is 



center to register the geocentric velocity frame with an accuracy of less than 0.5 km s in 
each spectrum. Translation from the geocentric frame to the Local Standard of Rest (LSR) 
is a simple calculation that depends on the direction and time of the observation. We use the 
same value for the Standard Solar Motion as that adopted by 21 cm observers, -|-20 km s~^ 
toward ccigoo = 18^*0, ^igoo = -|-30?0. 



Wc observed "standard" nebular sources at least every 90 minutes during survey obser- 
vations. All calibrations are tied to an absolute intensity measurement of the North American 
Nebula (NAN; NGC 7000) within a 50' beam at 850 R by Scherb (1981). Observations of 
NAN with WHAM yield the conversion: 1 R (km s~^)~^ = 684 ADUs for a 30 s exposure, 
for an estimated 800 R within our one-degree beam centered at ^ = 85?6, b — — 0?72. In the 
absence of extinction, 1 R corresponds to an emission measure (EM= J n1 dl) of 2.25 cm~^ 



pc for gas at a temperature of 8,000 K, a typical value for the WIM (Haffner et al. 1998, 
1999). 



By analyzing the nebular observations taken throughout each night at a variety of 
airmasses, we are able to compute the zenith transmission of the atmosphere at Ha for most 
nights of the survey. For nights where insufficient nebular observations exist for a specific 
determination, the atmospheric transmission is interpolated from an empirical fit to nights 
with well-determined transmission characteristics. We find that over the duration of the 
survey, the zenith transmission is described well by: 



where d is the day of the year, and A are the phase and amplitude of the seasonal trans- 
mission variance, and Tq is the yearly average transmission. The values for these parameters 




3.5. Intensity Calibration 



r = ro + Acos(— (ci-0)). 



(2) 
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derived from our calibration observations are To = 0.924, A — 0.026, and = 0.12. However, 
the root-mean-squared scatter of the real observations about this fit is large, 0.06, due to 
daily variations. With a transmission estimate available for each night of data, all intensities 
have been corrected to their intrinsic values above the atmosphere with the standard factor, 
rp-i/cosz^ where T is the zenith transmission of the given night and z is the zenith distance 
of an observation. Using Tq, the average correction increases from 8% near the zenith to 
17% near z = 60°. 

Two other small systematic intensity corrections have been applied. The first accounts 
for a small systematic distribution in the intensities of fully atmospheric-corrected nebular 
observations that appears to track the difference in tune pressures used for a set of obser- 
vations. We believe this arose due to the slight difference in tuning the elaton system when 
using a narrow source (the geocoronal Ha line) versus a broader- line source (diffuse nebulae) . 
This correction can increase or decrease observed intensities by up to 5%. The last correc- 
tion tracks the general degradation of the optical transmission of the system (typically due 
to dust on optical surfaces). Over the 24 months of the survey, this rate has been roughly 
constant at about 0.0103% per day, resulting in a gradually increasing correction of up to 
about 7.3% over the time the survey data were collected. 

Gaustad et al. (2001) find that our calibrated data match their independent determina- 
tion to better than 10%. 

3.6. Zodiacal Absorption Correction 

As a testament to the sensitivity of WHAM, one subtle and unexpected feature remains 
in many of the survey spectra after all atmospheric emission lines are removed. Figure 5 

shows a pair of evening "block-averaged" spectra where the Galactic emission is simple (i.e. 
well-described by a single component) and faint (Jhq; ~ 0.5 R). A representative fit to each, 
including the background. Galactic, and atmospheric components (§3.3), is also shown, as 
well as the residuals from these fits. A subtle negative excursion remains in the residuals of 
Figure 5b slightly blueward of the geocoronal line. This residual feature is evident at some 
level in 10-20% of the NSS blocks. Several pieces of evidence reveal that this feature is the 
solar Fraunhofer Ha absorption line scattered off of zodiacal dust. Most importantly, those 
blocks observed with small angular distances to the Sun projected along the ecliptic plane 
(quantitatively known as the solar elongation, e = A — A©) and closest to the ecliptic plane 
delineate the anomalous features in our maps. Such close correlation with this coordinate 
frame points strongly to the feature arising from zodiacal light, whose intensity is also best 
described in this coordinate system. Also, the location of the feature relative to the geocorona 



- 15 - 



(a fixed geocentric reference point near Vgeo — km s~^) changes from evening to morning 
observations. In evening observations the feature appears blueward of the geocorona, while 
in the morning it appears redward. This change in the location of the feature from evening to 
morning is consistent with the expected velocity shift of the solar spectrum due to the orbital 
motion of the reflecting dust along these two different lines of sight where e changes sign. 
The presence of this feature succinctly explains the deficiency of emission in evening spectra 
revealed by uncorrected channel maps centered at vlsr = ~50 km and in morning spectra 
by the i'lsr = +50 km s"^. However, it also suggests that these contaminated spectra should 
be wholly corrected to recover emission removed by the presence of this broad absorption 
feature over the whole spectral window. 

Because the magnitude of this effect is small (typically < 0.3 R) and only a small portion 
of the absorption feature remains uncontaminated by emission features, we are unable to 
reconstruct the absorption profile from the data alone. Instead, we adopt a simple, global 
approach to minimizing this effect in the final data set. First, we assume that the zodiacal 
absorption line has the same shape as the solar Fraunhofer line. An electronic version of the 
spectrum of the solar Ho; line from Delbouille et al. (1973) is used to construct a template 
absorption profile. Second, we assume that the depth of the absorption line scales linearly 
with the visual surface brightness of the zodiacal light. Using the data provided in Cox 
(2000), we derive a smooth relationship for the V-band surface brightness as a function of 
e. Combining this function with a empirically determined constant scales the absorption 
profile to best fit our data. Finally, we use the velocity model of Fried (1978) to estimate 
the velocity shift we should apply to our model absorption line as a function of e. We note 
that, prior to WHAM (Reynolds et al. 2003, in preparation), determinations of this final 
parameter have been fraught with large uncertainties, and studies have not agreed on the 
shape of the zodiacal rotation curve (Fried 1978; East & Reay 1984; Robley et al. 1985; 
Lcvasscur-Rcgourd et al. 2001). However, given the confusion of the gcocoronal emission 
and zodiacal absorption in the data under consideration, we believe that additional accuracy 
of these velocity shifts will not substantially aid the removal of this feature from the survey 
data. 

Even with this model, accurately removing this feature from the spectra is not com- 
pletely satisfactory. As Figure 5 shows, the geocorona and Galactic emission obscure most 
evidence of an absorption feature. As a result, the strongest residual of the zodiacal feature 
in the most heavily contaminated spectra only occurs near one wing of the geocoronal line — 
that on the opposite side of the Galactic emission. With only this small area significantly 
constraining our overall scaling of the model, we find that a single factor is insufficient to 
adequately correct for the absorption line. The final scale factor for the evening spectra 
(e < 0) is 2.5 times greater than for morning spectra. The origin of this asymmetry is still 
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unknown. The average strength of the hne added back in to the spectra is 0.2 R, with a 
maximum hne strength of 0.4 R for a few blocks taken particularly close to the Sun. Even 
with this correction applied to our final data set, special care should be taken when explor- 
ing survey data within 15 km s^^ of the geocentric rest wavelength of Ha because of the 
uncertainties associated with removing this zodiacal dust absorption feature. Of course, the 
removal of the much larger geocoronal emission feature also contributes to the uncertainty 
of the interstellar Ha spectrum near the earth's rest velocity. 

3.7. Maps 

Images (i.e. "maps") created from the survey data and presented here arc constructed 
by first integrating each original spectrum over the velocity band of interest. Then, using 
the coordinates for the observations and these integrated intensities, our irregularly gridded 
data is translated onto a regular grid by Delaunay triangulation (IDL routines TRIANGU- 
LATE and TRIGRID). For the data presented here, observations contaminated by bright 
stars (my < 5) are excluded from the interpolation and replaced with an average of uncon- 
taminated observations within one degree of the original pointing. 

4. SURVEY RESULTS 

There are a plethora of new discoveries throughout the interstellar medium that are 
arising from this Ha survey and from other emission lines being probed by WHAM. Some of 
the early results have already been pubhshed elsewhere (Reynolds et al. 1998; Haffner et al. 
1998; Tufte et al. 1998; Haffner et al. 1999; Reynolds et al. 1999; Lagache et al. 2000; Heiles 
ct al. 2000; Reynolds et al. 2001b; Haffner ct al. 2001; Reynolds et al. 2001a; Hansen et al. 
2002b; Tufte ct al. 2002; Hansen et al. 2002a). In this paper, we highhght the global results 
of the survey. 

4.1. Total Intensity Map 

Figure 6 shows a series of three maps (centered at different projection centers) of the Ha 
emission from the northern sky, integrated between Ulsr ~ 80 and +80 km s~^. In addition 
to the bright knots of emission from discrete H II regions in the plane, the general diffuse 
glow of the WIM is easily seen stretching to the Galactic poles. This material is rich in 
structure, with brighter filaments, loops, and regions superimposed on a fainter, smoother 
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(or unresolved) background. As §4.4 below discusses, the total emission is roughly described 
by a plane-parallel layer about the Galactic midplane. 

Due to the large dynamic range in the data, structure in the polar regions is lost in 
these all sky representations. Figure 7 focuses on only these regions, revealing that the 
diffuse Ha emission has detectable emission (and structure) everywhere in the sky. As more 
fully detailed in §4.4 below, very few beams in the survey have total intensities below a few 
tenths of a Raylcigh. A few directions with very low (or no) emission detected in these short 
30-second exposures have been studied in more depth by Hansen et al. (2002b). 

Table 2 gives a more quantitative representation of these integrated intensity Figures. 

4.2. Channel Maps 

The power of the WHAM-NSS is revealed in Figure 8. Here we present individual 
channel maps produced by using smaller integration ranges (Ai>lsr = 20 km s~^). The 
Galactic plane now breaks up into emission produced by individual spiral arms. In particular, 

the outer Perseus Arm is prominent at negative velocities (wlsr < ~30 km s~^) for i = 90° 
to 150° and at positive velocities (vlsr > +30 km s^^) for i = 210° to 240°. The inner 
Sagittarius- Carina Arm is present primarily at (vlsr > +30 km s~^) for £ = 0° to 50°. 

Many discrete features seen in H 1, especially those at intermediate velocities, are also 
seen in the WHAM-NSS (see §5). A few prominent structures, such as the 50°-long filament 
discussed in Haffner et al. (1998) rising near £ — 225° do not appear to have obvious neutral 
phase analogs. 

The spectral data from the WHAM-NSS is available through the Catalogues service at 
CDS (http: / / cdsweb.u-strasbg.fr/) . 

4.3. Longitude- Velocity Diagrams 

We present a few representative £-v diagrams in Figure 9. Although they highhght the 
major features of the Galactic velocity distribution of the ionized gas including spiral arm 
delineation and velocity extent, they do not show detailed structure at the level found in 
displays of H I and CO emission. This is primarily due to the lack of Ha components with 
narrow widths. The thermal broadening of the Ha line in gas with temperatures near 10*^ 
K combined with any non-thermal motions result in line widths typically greater than 25 
km s^^. 
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In addition, the effect of extinction by interstellar dust must be taken into account when 
using such diagrams to interpret global Galactic ionized gas distributions, especially within 
about 10° of the Galactic plane. Above 10°, these effects are small on average over the 
sky. However, very near the plane, some of the path lengths are known to be truncated by 
extinction, leading to potential confusion when interpreting a longitude-velocity diagram. 
We omit this particular latitude slice from Figure 9 as a result. 

4.4. General Quantitative Analyses 

Inspired by the discussion of the global properties of H I emission in Dickey & Lockman 
(1990), we present a similar analysis here for the Ha emission over the 3/4 of the sky covered 
by the WHAM-NSS. Figure 10 presents the distribution of the H« intensities as sampled 
by our one-degree beam over the sky. With the influence of discrete H II regions, the Ha 
emission spans over five orders of magnitude even before any corrections are made for dust 
extinction and dilution by our one-degree beam. The few brightest pointings (in particular 
near the Orion Nebula, M42, and H II regions in the plane toward the inner galaxy such as 
M20) measure over 1000 R. They are certainly not resolved by WHAM and have intrinsic 
Ha intensities much higher. On the faint end of the Ha intensity distribution, there is 
a very sharp cutoff between log = —0.75 and —0.35 (0.2 and 0.4 R). This transition is 
comfortably above the sensitivity of the 30-second integrations used for the survey (3(7 ~ 0.15 
R). As with the neutral gas, the minimum emission in Ha is not at the Galactic pole, but 
in the direction of the Lockman Window near £ = 152°, b = +52° (Jahoda et al. 1990; 
Lockman et al. 1986). The regions of lowest ionized and neutral gas emission are not exactly 
spatially coincident but are within a few degrees of one another. A one-degree beam within 
the Lockman Window and a region of even lower Ha emission nearby have been examined 
in detail by Hansen et al. (2002b). 

Figures 11 and 12 examine the distribution of total lua as a function of Galactic coordi- 
nates. The top panel of both figures plot each pointing in the survey while the bottom panel 
shows the median and spread (using the average deviation) within small bins (~ 1°) along 
the coordinate axis. Although using the median statistic does remove some of the influ- 
ence of anomalously bright regions, the imprint of a few very large emission regions persists 
nonetheless. In particular, the Orion-Eridanus superbubble covering the large region from 
i = 180° through 220° and 6 = 0° through b = -50° (Hcilcs et al. 1999) causes a substantial 
enhancement in the trend of the median and size of the spread in both figures. 

The sharp cutoff at low intensities in Figure 10 and the distinctive shape of the latitude 
profile of the Ha emission in Figure 11 can be explained to some extent by adopting a simple 
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model of the ionized gas as a uniform plane-parallel layer. In this case, the minimum path 
through the layer at the Galactic poles provides a well-defined minimum intensity for any 
direction while the increase of the emission at low b is mostly due to the geometrical effect 
of longer pathlengths through the layer. In such an idealized model, this cosecant effect is 
removed by plotting /HQsin|6| versus sin|6|, as shown in Figure 13. The top panel again 
plots every point in the WHAM-NSS and the bottom panel plots the median and average 
deviation, with the trend for the northern Galactic hemisphere {b > 0°) split from that of the 
south (6 < 0°). This split is useful for two reasons. First, the longitude coverage of the north 
is much more complete. Second, the Orion- Eridanus superbubble (seen as the enhancement 
between sin |6| = 0.1 and 0.4) only contaminates the southern data. 



5. DISCUSSION 

A full exploration of the rich detail in the WIM uncovered by the WHAM-NSS is 
beyond a single work. Here we comment on a few of the global observations highlighted by 
the presentation of the data in the various forms of §4. 

As also discovered by the recent deep Ha imaging surveys (Gaustad et al. 2001; Dennison 
et al. 1998), the distribution of the ionized gas in the Milky Way is characterized by both 
a pervasive background of faint emission and rich filamentary structure. Now, with the 
absolute calibration provided by WHAM and Galactic emission profiles free from atmospheric 
and solar system contamination, we can quantify the background accurately and, in many 
cases, isolate the brighter features spectrally. Haffner et al. (1998) presented one of the more 
striking features evident in Figures 6 and 8, a > 50° long filament stretching out of the 
Galactic plane near i = 225°, showing changes in its radial velocity along its length. Unlike 
most other distinct features in the ionized gas, this filament appears to have no counterpart 
in the neutral phase. Reynolds et al. (2001a) explore a large arc of emission above the Perseus 
arm near i — 135° (most visible in the lowest two velocity channels of Figure 8). When the 
NSS was completed, a similar arc was discovered south of the plane, directly opposite of the 
northern extension (Madsen et al. 2002). These massive structures were shaped out of the 
interstellar medium some time ago, but are now actively ionized by younger OB associations 
below them in the Galactic plane. Many of the intermediate-velocity complexes cataloged 
by Wakker (2001) also have associated Ha emission revealed in the WHAM-NSS. Haffner 
et al. (1999) discuss the character of one of the largest revealed in our survey, Complex K. 
The discovery of this associated ionized gas will not only aid the characterization of the 
intermediate- and high- velocity clouds (Wakker et al. 1999), but also provides one of the 
only direct probes of the ionizing flux in the halo (Tufte et al. 1998; Bland-Hawthorn & 
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Maloney 1999; Tufte et al. 2002). Finally, the survey has revealed several new H II regions 
around early B and hot subdwarf stars (Haffner 2001). These discoveries could be important 
by helping to characterize a population of sources capable of providing ionizing flux away 
from the heavily absorptive Galactic disk. 

The kinematics of the WIM are broadly summarized by Figure 9. Two general trends 
arc worthy of special note. First, the kinematic imprint of spiral arms in both the inner 
and outer Galaxy is clearly seen even 30° off the plane. Above this height, only a few 
discrete features that stretch to large heights above arms and intermediate- velocity complexes 
still show substantial non-local velocity patterns. The second general trend is that at high 
latitudes the Ha emitting gas is clearly biased toward negative velocities. The emission at 
\b\ = 60° on both sides of the Galactic plane shows an asymmetric distribution about the 
solid hue that denotes i'lsr = 0. This trend is also seen in the neutral gas toward the 
Galactic poles (Blitz et al. 1999). Figure 14 shows the spatial distribution of this effect by 
presenting maps of the average intensity- weighted velocity (J I^avdv/ J luadv) over both 
Galactic poles. In particular, it highlights the fact that deviations from zero velocity occur 
over coherent regions of rising and falling gas. 

Even when the localized impact of bright H II regions is minimized (§4.4), a plane- 
parallel model turns out to be only approximate in describing the global WIM, as is also 
the case for the neutral gas (Dickey & Lockman 1990). The dotted line in Figure 13 delin- 
eates the overall median value of log/HQsin|6| for the survey, —0.17, corresponding to an 
/HaSin|6| = 0.68 R. Extinction by dust likely creates the observed decrease in /HQsin|6| at 
sin |6| < 0.05 (i.e. |6| < 3°). Aside from this expected effect, instead of being constant, the 
Galactic emission decreases by a factor of two in the northern hemisphere from sin |6| = 0.2 
through 1.0. Even with the limited longitude coverage in the south, the trend appears to 
be similar. As pointed out by Dickey & Lockman (1990) to explain the similar effect in the 
H I distribution, the presence of the Local Bubble (Maiz-Apellaniz 2001; Smith & Cox 2001; 
Cox & Reynolds 1987), a localized region of very low density ionized gas around the Sun, 
likely causes at least some of the decrease seen here. However, the much larger scale height 
of the Ha emitting gas (Haffner et al. 1998; Reynolds 1997) should reduce the imprint of this 
local feature compared to its importance in the same analysis of the neutral gas. Two other 
effects unique to the ionized gas emission may also contribute to the slope seen here. As dis- 
cussed in Haffner et al. (1998) and Reynolds (1997), the vertical electron distribution derived 
from Ha observations appears to be fairly well described by a two-parameter exponential. 
However, these parameters may be determined primarily by local star formation so that, in 
particular, massive spiral arms may have distinctly different sets of parameters compared 
to the local neighborhood in which the Sun resides. In this scenario, low-6 sightlines will 
be preferentially elevated, resulting in a decrease in /nasin \b\ with increasing sin \b\ as the 



- 21 - 



local parameters become increasingly influential in the total emission along the line of sight. 
Finally, the increase at lower latitudes may be due in part to contamination by Ha scattered 
by interstellar dust that is far from the emitting source (Wood & Reynolds 1999). 

6. SUMMARY 

We have completed the first velocity-resolved survey of Ha emission over the northern 
sky. The WHAM-NSS sheds new hght on the distribution and kinematics of the WIM, 
allowing comparative studies on a global scale with other phases of the interstellar medium. 
The survey reveals that: 

• Ionized gas in the Milky Way is ubiquitous, with a pervasive yet faint background 
{Ina > 0.1 R) detectable in nearly every direction in the northern sky. 

• The total Ha emission is roughly described by a plane-parallel layer about the Galactic 
midplane having /naSin |6| = 0.68 R, although there are noticeable systematic devia- 
tions from this average value. 

• The vertical distribution of the ionized gas is more extended than that of the neutral 
gas (Haffner et al. 1998; Reynolds 1997). 

• The average kinematics of the ionized gas follow the velocity patterns of the nearby 
spiral arms and the imprint of the arms is able to be traced to a height of |6| = 30° in 
most cases. 

• At high latitudes, the kinematics of the ionized gas is preferentially biased toward 
negative velocities. 

• Many of the anomalous intermediate- velocity features (|vlsr| < 100 km s~^) discovered 
in H I surveys also contain detectable amounts of ionized gas (Haffner et al. 1999). 
WHAM's detection of Ha emission from high-velocity clouds (Tufte et al. 1998) is not 
covered by this survey. 

The WHAM-NSS along with the new deep Ha imaging surveys by (Gaustad et al. 2001) 
and (Dennison et al. 1998) open a host of new opportunities to study the WIM. WHAM is 
continuing to explore the ionized gas of the Milky Way by extending observations to higher 
velocities and adding additional diagnostic emission lines in regions of interest. 

The full WHAM-NSS is available at http://www.astro.wisc.edu/wham/. 
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Table 1. Faint Atmospheric Lines near Ha 





^geo 


FWHM 


Relative 


Line 


(km s~^) 


(km s~^) 


Intensity 


1 


-148.5 


10 


0.54 


2 


-131.6 


10 


0.24 


3 


-88.6 


10 


0.20 


4 


-61.7 


15 


0.33 


5 


-42.0 


10 


1.00 


6 


-25.1 


10 


0.30 


7 


+32.4 


10 


0.30 


8 


+40.7 


10 


0.65 


9 


+53.9 


10 


0.39 


10 


+72.6 


15 


1.60 


11 


+95.9 


10 


0.36 


12 


+120.6 


10 


0.29 
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Table 2. Integrated Ho; Emission between t'LSR= —80 and +80 km s 



Position 


-'Ha 


o 




xJl igjiiiai 




h 










76? 189 


-46?682 


1.204 


0.000 


1000 


1.204 


77?619 


-46?682 


1.285 


0.034 


1000 


0.000 


79?049 


-46?682 


1.047 


0.035 


1000 


0.000 


80?479 


-46?682 


0.700 


0.034 


1000 


0.000 


81?569 


-45?832 


0.668 


0.034 


1000 


0.000 


80?169 


-45?832 


0.634 


0.037 


1000 


0.000 


78? 759 


-45?832 


0.867 


0.035 


1000 


0.000 


77?349 


-45?832 


1.205 


0.035 


1000 


0.000 


75?949 


-45?832 


0.993 


0.035 


1000 


0.000 


75?309 


-44?981 


1.108 


0.034 


1000 


0.000 



^1 R = 10^/471 photons cm-^ s^^ sr-^ = 2.4 x 10"^ 
ergs cm~^ s~^ sr~^ at Ha. 

^Formal propagation only. This error does not re- 
flect systematic issues discussed in the text. If hsted 
as zero, the original observation is likely contaminated 
by a bright star (see text). The intensity listed in the 
third column (/Ha) is then an average of the nearest 
neighbors within one degree of the original. 

^The numeric block designation for this observation. 
Blocks of data arc taken sequentially in time and arc 
physically close on the sky (see text). Systematics re- 
maining in the data are thus often linked to this pa- 
rameter. 

^If the third column represents the average of nearest 
neighbors (i.e. the fourth column is zero; see note h 
above), this column contains the original (unmodified) 
intensity for reference. 

Note. — The complete version of this table is in the 
electronic edition of the Journal. The printed edition 
contains only a sample. 
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Fig. 1.— The WHAM instrument. The top photo shows WHAM at Kitt Peak National 
Observatory, AZ. A schematic of the instrument is displayed in the bottom panel. The 
dashed line denotes the path of an on- axis ray through the system. As shown in the photo, 
both flat mirrors in the siderostat rotate independently, allowing all-sky pointing. 
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Fig. 2. — Etalon transmission functions. These curves depict the theoretical fraction of 
hght transmitted through the (a) high-resolution and (b) low-resolution Fabry-Pcrot ctalons 
used in WHAM. This example shows the transmission functions for a single incident angle 
with the etalons "tuned" (§2.1) near Ha. When used in series, this pair of etalons produces 
the composition transmission curve shown in (c). A theoretical filter profile with FWHM 
fa 20 A is shown as the dashed fine, (d) plots the logarithm of the composite etalon and 
filter transmission function, showing the detailed structure of the residuals from the order 
suppression. 
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Fig. 3. — Survey beam sample. Each circle represents a one-degree pointing of the WHAM 
survey. Pointings of the same shade of gray belong to the same observational "block" (see 
§2.2). 
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Fig. 4. — WHAM data reduction. An unprocessed ring image is depicted in (a). This image 
is "ring-summed" (see §3.1) to produce (b). The fully- processed, calibrated spectrum in (c) 
contains only Galactic Ha emission toward i — 124?5, b — 1?7. 
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Fig. 5. — Zodiacal Absorption Feature. Data and reduction fits are shown for two "block- 
averaged" spectra. Uncalibrated CCD counts are plotted logarithmically against a velocity 
frame with an arbitrary, but standardized zero-point, (a) is the average spectra of the 41 
pointings in the block centered near i — 167?5, b — 32?2. (b) is the average spectra of the 42 
pointings in the adjacent block centered near i — 174?4, b — 32?2. The data are plotted as the 
-I- symbols. Each plot shows the total fit (solid line) as well as the individual components 
(labeled above) used to remove the various contaminating elements discussed in §3. The 
residuals of each fit are plotted in the bottom panel of each figure. The increased scatter 
in the residuals around v = 30 km s~^ is expected statistically due to the increased signal 
of the geocorona, whose peak contains more than 20 times more counts than other portions 
of the spectrum, (b ) reveals a systematic residual present in many blocks near the ecliptic 
plane. The feature is undetectable in (a) because it is 5° further from the ecliptic plane. §3.6 
details our attempt to minimize the effect of this hidden feature, a portion of the zodiacal 
Praunhofer absorption line. 
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Fig. 6. — WHAM-NSS: Total Intensity. The integrated Ha intensity between wlsr = —80 
and +80 km s~^ is mapped as a function of Galactic coordinates. These three Hammer- 
Aitoff projections are centered at 6 = 0° and £ = 120°, 0°, and 240°, respectively from top 
to bottom. Dashed lines are spaced 30° apart in longitude and 15° apart in latitude. 
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Fig. 7. — WHAM-NSS: Total Intensity (Polar View). These stereographic projections map 
the total integrated Ha intensity between wlsr = ~80 and +80 km s~^ for the (a) northern 
(6 = +40° to +90°) and (b) southern {b = -40° to -90°) Galactic polar caps. Green dashed 
lines are separated by 60° in longitude and 10° in latitude. The grid line for i = 180° is 
at the twelve-o'clock position. The solid line marking the border of data in the southern 
projection denotes S = -30°, the hmit of the WHAM-NSS. 
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Fig. 8. — WHAM-NSS: Channel Maps. This series of maps shows the kinematic distribution 
of Ha emission in the WHAM-NSS. Each map is integrated over the 20 km s~^ range hsted 
at the top of the map. Each map is scaled independently through histogram equalization to 
highlight the substructure within the channel. 
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Fig. 9. — WHAM-NSS: Longitude- Velocity Diagrams. Each map displays the average 
longitude-velocity diagram for a 5-degree wide strip around the latitude listed above the 
map in steps of 10° of latitude. Each map is centered around i = 120°. The left column 
plots positive latitudes {b = +10° through -|-60°) and the right plots negative latitudes 
[b = —10° through —60°). Both increase in distance from the Galactic plane toward the 
top of the figure. Since 6 = 0° is heavily influenced by dust extinction, it is not shown (see 
text). The intensity of a point in the diagram is represented by logarithmically scaling the 
spectral values to a full-color table. For rough comparison to integrated intensities, a single 
component with FWHM = 30 km s~^ and a peak value of 2.7 x 10~^ R (km s~^)~^ has a 
total Ha intensity of 1 R. 
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Fig. 10. — Distribution of Intensities. This histogram plots the distribution of the log/na 
values in the WHAM-NSS. Notice that the steep drop-off near log lHa= —0.5 (~ 0.3 R) is 
above the sensitivity of our instrument (3(7 ~ 0.15 R, denoted by the dotted line). Within a 
one-degree beam, most directions have intensities within a decade of 1 R (EM ~ 2 pc cm~^). 
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Fig. 11. — Intensity versus Galactic Latitude. The top panel plots the logarithm of the 
total intensity of every point in the WHAM-NSS as a function of Galactic latitude. Since 
the survey is regularly gridded in latitude, a small random perturbation (up to ±0.85°) is 
applied to the true latitude of each point in the plot to avoid striation. The bottom panel 
plots the median value within each of the 0.85° latitude strips from the survey. The vertical 
extent of each plotted bar is determined from the average deviation about the median within 
that bin. The larger deviations at 6 = —10° through —50° are due to the presence of the 
Orion-Eridanus superbubble (a brighter than average local feature) at these latitudes. 
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Fig. 12. — Intensity versus Galactic Longitude. The top panel plots the logarithm of the 
total intensity of every point in the WHAM-NSS as a function of Galactic longitude. The 
bottom panel plots the median value within one-degree longitude strips from the survey. The 
vertical extent of each plotted bar is determined from the average deviation about the median 
within that bin. These median values are only shown for the regions of the WHAM-NSS 
that are complete in latitude. 
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Fig. 13. — /HaSin|6| vs sin|6|. If Ha emission were to arise from a plane-parallel slab, 
this representation would show intensity constant at all latitudes. The two tracks in the 
bottom panel divide the data into positive (black) and negative (gray) Galactic latitudes. 
Emission near sin |6| = is likely influenced by extinction. The dotted line denotes the 
median Jho sin |6| value for the survey pointings. 
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Fig. 14. — Velocities at high-latitudes. These maps mirror the projection aspects of Figure 7 
but display intensity- weighted velocity averages (see §5), which highlight a negative velocity 
preference for ionized emission toward both Galactic poles. 



